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ABSTRACT 

Motivated by detections of hypervelocity stars that may originate from the Galactic Center, we 
revist the problem of a binary disruption by a passage near a much more massive point mass. The 
six order of magnitude mass ratio between the Galactic Center black hole and the binary stars allows 
us to formulate the problem in the restricted parabolic three-body approximation. In this framework, 
results can be simply rescaled in terms of binary masses, its initial separation and binary-to-black hole 
mass ratio. Consequently, an advantage over the full three-body calculation is that a much smaller 
set of simulations is needed to explore the relevant parameter space. Contrary to previous claims, 
we show that, upon binary disruption, the lighter star does not remain preferentially bound to the 
black hole. In fact, it is ejected exactly in 50% of the cases. Nonetheless, lighter objects have higher 
ejection velocities, since the energy distribution is independent of mass. Focusing on the planar case, 
we provide the probability distributions for disruption of circular binaries and for the ejection energy. 
We show that even binaries that penetrate deeply into the tidal sphere of the black hole are not 
doomed to disruption, but survive in 20% of the cases. Nor do these deep encounters produce the 
highest ejection energies, which are instead obtained for binaries arriving to 0.1 — 0.5 of the tidal radius 
in a prograde orbit. Interestingly, such deep-reaching binaries separate widely after penetrating the 
tidal radius, but always approach each other again on their way out from the black hole. Finally, our 
analytic method allows us to account for a finite size of the stars and recast the ejection energy in 
terms of a minimal possible separation. We find that, for a given minimal separation, the ejection 
energy is relatively insensitive to the initial binary separation. 

Subject headings: Galaxy: Center, Galaxy: halo, Galaxy: kinematics and dynamics, Galaxy: stellar 
content, Binaries: general 



1. INTRODUCTION 

Hypervelocity Stars (HVSs) are stars with a velocity 
exceeding the escape velocity of the Galaxy. Currently, 
16 such stars have been observed, 15 of whic h are thought 
to o r iginate from our Galactic Centre ([Brown et al.l 
l200l 120061 . 120071 . [2001 IHirsch et al.ll2005l) and one from 
the Large Magellanic Cloud (|Edelmann et alJ l2005h . 
All have been observed with radial velocities between 
300 km/s and 800 km/s and almost all are located over 
50 kpc away. For stars whose likely origin is within the 
Galaxy and taking the galactic potential into account, 
this translates to velocities of over 1000 km/s from the 
bulge. The observational strategy is such that most of 
the discovered HVSs are faint B stars. They have es- 
cape velocities from the surface of the order of 600 km/s, 
well below the ejection velocity of 1000 km/s. Thus, the 
standard mechanisms for producing high velocity run- 
away stars, such as star-scattering and explosion as a 
supernova of one component of a binary cannot work. A 
dynamical interaction with a massive compact object is 
likely involved. 

In this paper, we adopt one of the leading models for 
the formation of HVSs: the breakup of a binary as it 
appro aches the black hole in the Galactic Center (|Hillsl 
1988). Simple analytical arguments can be made to show 



the potential of this model to explain HVSs. 

If the binary of total mass m has separation a, then 
tidal forces from the black hole overcome the binary's 
mutual gravitational forces at the tidal radius r t — 
a(M/m) 1 / 3 , where M is the mass of the black hole. The 
relative velocity of the binary components is of order 
vq = (Gm/a) 1 / 2 . If the binary approached the black hole 
with negligible energy, its center of mass moves at the 
tidal radius with velocity of order vbh = (GM/rt) 1 ^ 2 — 
voiM/m) 1 / 3 relative to the black hole. There are three 
ways to estimate the energies of the individual compo- 
nents of the binary, assuming that they arrived with neg- 
ligible total energy. It is instructive to consider all three: 

• Kinetic energy: Adding or subtracting the rela- 
tive velocity of the components, vq, to the velocity 
around the black hole, results in an additional ki- 
netic energy of order «o u bh ~ v 2 (M /m) 1 / 3 . 

• Gravitational potential energy: The displace- 
ment of order o in the position of each component 
of the binary, at a distance of about r t from the 
black hole, results in a change in gravitational en- 
ergy of GMa/r 2 ~ ^(M/m) 1 / 3 . 

• Work: The energy of each of the component in the 
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black hole frame, is changing only due to mutual 
forces between the binary components. The force 
is of order Gin /a 2 and the length, in the black hole 
frame, over which it acts is tv Therefore the work 
is Gmr t /a 2 ~ ^(M/m) 1 / 3 . 

All these estimates provide an energy of order 
vliM/m) 1 ' 3 . If the binary dissolves, one component of 
the binary stays bound to the black hole, the other es- 
capes, with velocities at infinity 1 of vq{M /m) 1 ^ . The en- 
counter with the black hole therefore allowed for a larger 
velocity by a factor of (M/m) 1 ^ 6 than the orbital veloc- 
ity of the binary. For the parameters of our Galaxy and 
stars, (M/m) 1 / 6 ~ 10, allowing ejections with thousands 
of kilometers per second. 

This is the theoretical framework in which iHillsl (|1988f ) 
predicted the existence o f HVS s. L ater, it was dis- 
cussed by lYu fc TremaTnel (|2003h and iGould fc Quilled 
(2003). After the observational discovery, many 
more papers on the subject appe a red ([Gualandris et al.l 
2001 IGinsburg fc Loebl [200 6a ,b: B romley et alj 120061: 



Sesana et al.l 120071; iPerets et all 120071 : iKenvon et all 
2008T ) . aiming to predict the properties of the 
ejected and/or captur ed stars. Ot her p apers, (e.g. 
IGinsburg fc Loebll2007l : lAntonini et"aTll2009l ). focused in- 
stead on the fate of binaries that are not dissolved, but 
that in fact coalesce. The investigations so far have used 
three-body simulations or analytic methods that relied 
on results from three-body simulations. Consequently, 
only a limited set of parameters (e.g. for the binary mass 
ratio) have been explored. 

In this paper, we show that this problem can be in- 
vestigated with methods related to those used in the 
study of asteroids in the Solar System, which exploit 
the enormous disparity in mass between the bodies in- 
volved. Specifically, we can formulate it in terms of a 
restricted 3 body problem, i.e. the motion of a single 
massless particle under the influence of external time de- 
pendent forces. Our treatment is valid as long as the 
binary components are closer to each other than to the 
black hole. Since r t = a(M/m) 1 / 3 our approximation re- 
quires (M/m) 1 / 3 1, which is fully satisfied in the case 
of the black hole at the Galactic Center and a binary of B 
stars, (M/m) 1 / 3 = 100. With the advantages of this an- 
alytic method, we can reach general conclusions that do 
not depend on the physical properties of the system, such 
as masses and binary semi-major axis. Moreover, the or- 
bit integration is faster and more stable, allowing us to 
handle more easily cases of close encounters between the 
bodies. 

In <J2]we outline the formulation of the 3-body problem 
in terms of a restricted parabolic problem. In $3] we use 
the restricted radial problem to describes binaries that 
penetrate deep into the tidal radius. The radial problem 
has a singularity at the time that the binary encounters 
the black hole, and we use the results of the parabolic 
problem to pass the singularity and continue to describe 
the evolution of the binary on its trajectory away from 
the black hole. In <2]we compare our results with 3-body 

1 We ignore the galactic potential in this paper. By "the velocity 
at infinity" we mean the velocity of the object once it escaped the 
gravitational potential of the black hole, but did not yet climb out 
of the potential of the rest of the galaxy. 



numerical integrations, and find excellent agreement. We 
then use the parabolic and radial formalism to investigate 
the probability of dissolving a circular, planar binary and 
to obtain quantitative estimates of the ejection velocities. 
We outline these results in $5l 

2. THE PARABOLIC RESTRICTED 3 BODY PROBLEM 

2.1. The orbit 

In the case of interest M/ra 3> 1 and the equation of 
motion for each of the binary members reads 
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where ri and r 2 are the respective distance from the 
black hole. We therefore obtain that the equation for 
the distance between the two f = r 2 — ri is 
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Now, we assume that the two masses are much closer to 
each other, and to some effective parabolic trajectory r m , 
than each of them to the central mass M. As we already 
noticed in our Introduction, this condition is actually 
enforced, up to the tidal radius, by the requirement of 
a large mass ratio M/m ^> 1. Deep inside the tidal 
sphere, this assumption may break and we will discuss 
the relevant conditions in i)3.21 

We define the effective trajectory, r m , to be the 
parabolic orbit set by the position and velocity of the 
center of mass of the binary when it is far away from the 
black hole. The trajectory r m is therefore described by 



2r„ 



1 + cos / 



(4) 



where r p is the distance of closest approach, and / the 
angle from the point of closest approach. The angle /, 
know as the true anomaly, is a function of time, but one 
has analytically only the time as a function of /: 



t = 



^\/5 tan(//2) ( 3 + tan2 (// 2 )) ■ ( g ) 

Note that the actual center of mass will not generally 
move on this orbit, since the total force on the binary 
does not equal the force that would act if the binary was 
a single body. 

Linearizing the first two terms of equation ([3]) around 
r m , we find that the zero orders cancel out, and we obtain 

GM GM Gra. 
r = =- r + 3— r-(rr m )r m ^-r. (6) 

m m 

Rescaling the distance between the bodies by 
(m/M) 1 / 3 r p and the time by yJr^/GM, we can re- write 
eq. ([5]) in terms of the dimensionless variables 2 r and t 

2 Later, for the radial problem O, we will rescale our variables 
differently, keeping the same symbols. Therefore, in the whole 
paper, r and t should be understood as dimensionless, with a nor- 
malization that should be clear from the context. 
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^) [-r + 3(rr m )r m ]-4- (7) 



Since f m = (cos/, sin/, 0), r m /r p — 2/(1 + cos/) and 
we set r = (x, y, z), explicit equations in terms of dimen- 
sionless Cartesian coordinates read: 

.. (1 + cos/) 3 

£ = [—x + 3{xcosj + ysmf) cos j\ 



(x 2 + y 2 + z 2 ) 3 / 2 



y = g [-y + 3{x cos / + y sin f) sin /J 

2/ 



(1 + cos/) 3 
z = z 



(x 2 + y 2 + z 2 fl 2 

z 



t = 
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(9) 
(10) 



(x 2 + y 2 + z 2 ) 3 / 2 ' 
tan(//2)(3 + tan 2 (//2)). (11) 



Equations (I51)- ([TU1) are the equivalent of the Hill equa- 
tions ()Hilll 118861 ). for the parabolic rather than circular 
problem. Since equation pip is implicit, for numerical 
applications it may be preferable to use its differential 
form 

/ = V2(l + cos/) 2 /4. (12) 

2.2. Free Solutions 

Just like the Hill equations, equations ©-(US]) have 
free solutions, those that ignore the interaction term 
r/r 3 . These can be found mathematically, but physical 
intuition facilitates a swift solution. Since we have a set 
of three linear differential equations of the second order, 
all solutions are a linear combination of six independent 
solutions. Each could be physically obtained by taking 
the difference between an orbit infinitesimally close to a 
parabolic orbit and the parabolic orbit itself. We list the 
solutions below, stating which orbital element has been 
varied. 



• Variation in the argument of periapsis 
sin / cos / 
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1 + cos / 



y 



1 + cos / 



(13) 



• Variation in the time of periapsis 

x = — sin/ , y= 1 + cos/. (14) 

• Variation in the periapsis distance 

x = 2 — cos/ , y — — cos/tan(//2). (15) 

• Variation in the eccentricity at fixed periapsis 

x = (8 + 12cos/)tan 4 (//2), (16) 
35 sin / - 2 sin(2/) + 3 sin(3/) 



y 



(1 + cos/)* 



• Rotation around the apsidal line 

2 sin/ 



x = y = 



z 



1 + cos / ' 



• Rotation around the latus rectum 

2 cos/ 



x = y = , z 



cos/ 



(17) 



(18) 



In the above expressions, / is a function of t as given 
by equation (jlip . The first four solutions are planar, i.e. 
z = 0, while the last two solutions are one dimensional, 
x = y = 0. It is easy to check the validity of these 
expressions by substituting them into equations (f5|)- (fTU)) . 
and using eq. (fT2|) . 

2.3. The Energy 

We are ultimately interested in the fate of a star in 
a binary, following its encounter with the black hole. If 
the binary - approaching the hole on a parabolic orbit- 
is torn apart, a star can either become bound to the 
black hole or be ejected from the system. To distinguish 
between these two possibilities, we calculate its energy as 
a function of time, including the negative gravitational 
energy due to the black hole. Initially, at large distances, 



the specific energy of one member is simply ■ 



After 



the binary disruption, the analytical arguments in our 
Introduction suggest that its energy is larger by a factor 
of (M/m) 1 / 3 3> 1. We thus neglect the term due to the 
self-gravity of the binary. In addition, for (M/m) » 1, 
the change in energy of the massive black hole can also 
be neglected and the total energy of mi reads 

E 1 ^-^l +mil r 1 f/2. (19) 

ri 

To zero order, we can replace ri by r m , but since r m is 
a parabolic orbit that zero order energy vanishes. The 
first order terms are 



GMmi . 
E i = ( r i 



mir m (ri 



or, using our rescaled variables, 
GMm\m,2 ( m\ 1 / 3 



E, = -- 



M 



-rrv 



(20) 



(21) 



Since in our limit the total energy of the system is zero, 
the energy of the other body, ttt.2, is simply E2 = —E\. 
In the following, it is useful to define the penetration fac- 
tor D = r p /r t . In terms of our dimensionless Cartesian 
coordinates eq (f2"Tj) is given by 



£1 



(1 + cos/) 5 



Grnxm^ 

— -Cj2 — 

aD 

(xcosf + ?/ sin/) - 



M\ 1/3 
m J 



(22) 



sin fx + (1 + cos f)y 
71 . 



If the binary dissolves, this energy tends to a constant, 
since the body is eventually moving only under the con- 
servative force of the black hole. Mathematically, this 
means that the first term on the right-hand side of eq. ([7]) 
now dominates, and the problem is linear. 
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The negative of a solution is thus a solution. But the 
energy is also linear in the coordinates. Therefore, a 
body starting with 7r phase difference will have the same 
final energy in absolute value but opposite in sign. This 
is independent of the mass of the star. The important 
consequence is that, of the disrupted cases, half would 
have the heavier object bound and the lighter escaped, 
while half would have the opposite. These findings are at 
odds with those of Bromley et al (2006) who find that, for 
large mi/ma ratios, the lighter object usually becomes 
bound. 

We also note that, to this lowest order, the z compo- 
nent of the binary motion has no effect on the energy, 
thus it does not determine whether a body would be 
ejected. 

Finally, we explicitly write the energy, when we can 
neglect the interaction term between the two stars. The 
solution in this "free" regime 3 , is a linear combination 
of equations (flU|) - (fT5)) . For each free solution the energy 
is a constant of motion. Therefore the total energy is 
a linear combination of these constants. However, any 
solution constructed from two bodies in infinitesimally 
close parabolic orbits has zero energy. Only eq. (Tl6|) , that 
describes the relative orbits of two bodies where one with 
e ^ 1, gives a finite constant energy when substituted 
into eq. 

E 1 = -E 2 = -^^(^-) 1/3 20 B, (23) 
a D \ m J 

where B is the coefficient of eq. (|16[) in the linear expan- 
sion. 

3. DEEP PENETRATORS AND THE RADIAL RESTRICTED 
3 BODY PROBLEM 

The previous section allowed for arbitrary penetration 
factors D. The presence of the factor 1/D in eq. (|23[) 
seems to suggest that for a given initial star separation 
a arbitrary large energies can be attained for binary tra- 
jectories penetrating far into the tidal disruption region, 
(D <C 1). However, as we show here, this is not the case. 
We proceed to investigate this case of deep penetrators 
by considering a different, even simpler, limit of the equa- 
tions. The periapsis distance of the center of mass, r p , 
becomes now irrelevant and instead the center of mass 
can be taken to move on a radial orbit given by 



/ 9GMt 2 
V 2 



1/3 



(24) 



We can now scale the distance to the initial semi-major 
axis of the binary, a, and the time to the inverse an- 
gular frequency of the binary with that semi-major axis 
\J Gmj a? . Taking the direction of the binary center of 
mass with respect to the BH to be f m , equation |6]) sim- 
plifies to 



2 



3(rf„ 



r 



(25) 



3 Even in the non free regime, where the gravitational forces 
between the two stars are important, one can expand the solutions 
in terms of the free solutions, except that the coefficients will be 
time dependent. Our expression of the energy will still be given by 
B, except that B now is time dependent. 



In Cartesian coordinate, where x is the direction of the 
binary's center of mass, we have: 



X m 2 * {x 2 + y 2 + z 2 f' 2 ' 



9t 2 ' 
2 



(x 2 + y 2 + z 2 ) 3 / 2 



' 9t 2 (x 2 +y 2 + z 2 ) 3 / 2 ' 
The energy, as given by equation (|20|) is 



(26) 
(27) 
(28) 



E 1 



-E, 



Gm Y m 2 fM\ 1/3 f2^ 2/Z 



|t|- 4/3 (x + 3ii). 
(29) 



3.1. Special Solutions to the radial restricted problem 

3.1.1. Homogeneous Collapse 
It is simple to verify that the following 

x(t)=±(^j 1/3 \t\ 2 / 3 , y(t) = z(t) = 0, (30) 

is a solution to our equations. This solution is the analog 
of a homogenous collapse. In these solutions, the small 
masses mi and m 2 accelerate towards each other, by a 
combination of tidal forces and mutual acceleration, at 
the same rate that their center of mass is accelerating 
towards the central mass. These solutions, just like the 
stationary solution at the L\ and L 2 points in the cir- 
cular restricted problem, and just like the homogeneous 
collapse of dust, are unstable for deviations in the x- 
direction. 

3.1.2. Free Solutions 

In the limit of large separation between the small 
masses or close to t = where the distance to the black 
hole is small, tides dominate over the mutual gravity, and 
the problem become separable in x, y and z. It admits 
the following solutions: 



x(t)^A x \t\-^ + B x \t\ 4 / 3 , 
y{t) = A y \t\^ + B y \t\ 2 '\ 
z(t) = A z \t\ 1 / 3 + B z \t\ 2 ^. 



(31) 



These are the equivalent of the free solutions given for 
the parabolic restricted problem in equations (fll{l) - (fT5)) . 
Just like those, they can be given physical interpreta- 
tions. The A x solution describes two particles that have 
the same trajectory, but are slightly separated in time. 
The By and B z solutions, describe particles going on 
slightly different radial paths, each of zero energy. The 
B x solution describes the relative orbits of two particle 
going on the same radial path, but with slightly different 
energies. The A y and A z describe the relative orbits of 
particles with slightly different angular momentum. 

Note, that both exponents of the y and z terms are 
positive, while one exponent of the x term is negative. 
This means that at times close to zero, the particles are 
very close in y and z but are separated in x: y(t — > — Q) = 
z(t -> -0) = while x(t -0) = ±oo. 
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As in the parabolic case, one of the coefficients of the 
solution is related to energy. Substituting the equation 
for the x component and its derivative into eq. (|29[) we 

get, 

E 1 ^^^( A -i) 1/ \(2/9)^B x . (32) 
a \ m J 

This can be regarded has the equivalent of cq. (|2"3"1) in 
the limit of D — ► 0. It is evident, that the energy in this 
limit is finite, and the divergence suggested by eq. (|23[) 
is not real. Moreover, as we will show later, the highest 
ejection velocities are not obtained in this limit. 

3.2. The t = singularity 

From the discussion above, we learned that, beside a 
set of zero measure of initial conditions, the binary com- 
ponents acquire an increasingly large separation as they 
approach the large mass M. It is therefore tempting 
to conclude that for deep penetrators, binaries always 
dissolve. In fact, the two components of the binary ap- 
proach each other again at later (positive) times. Thus, 
to know the final outcome of the scattering, we should 
overcome the singularity and follow the orbits of the light 
bodies beyond the time of the encounter with the black 
hole. 

The difficulty arises from the assumption of a purely ra- 
dial orbit for the binary center of mass. A deep (D -C 1) 
parabolic orbit would parallel closely the radial one and 
get around the black hole smoothly. Once across the 
singular region, the radial equations are valid and the 
integration of the orbit can be resumed. We further no- 
tice that when the binary is well within the tidal ra- 
dius, \t\ <C 1 the mutual gravity of its members can 
be neglected and their separation follows a free solu- 
tion. Therefore, the resolution of our problem ultimately 
lies in finding which free parabolic orbit reduces -in the 
deep penetration approximation- to the free radial or- 
bit that an approaching binary is following, while still 
farther away from its periapsis, \t\ ^> D 3 / 2 . 

This translates in expanding eqs. (|13|) - (fl8l) around 
/ = — 7r and comparing them with eqs. (|31|) . Interest- 
ingly, we find that there is a one to one correspondence: 
solutions with coefficients A x , A y and A z correspond to 
eq. (fT4|) . eq. (fl"5f and eq. (|T7|) respectively. Likewise, 
the solutions with B x , B y and B z correspond to eq (|T6|) . 
cq. (fT3|) and eq. (fl8jl respectively. This is as expected 
given their physical interpretation. To make an exam- 
ple, we expand eq. (|11[) and the x component of eq. (| 16[) 
multiplied by some coefficient B around / = — n. We 
obtain 

x = -6 4 / 3 \t\ 4 / 3 B. (33) 

This is the behavior of the radial solution with coeffi- 
cient B x . In fact, comparing the two solutions for x and 
accounting for the different spatial and temporal dimen- 
sions, we get 

B x = -6 i/3 D- 1 B. (34) 

With this relation, the two equations (f23f and |32|) are 
equivalent as we expect. 

From this comparison, we also learn the behaviour of 
the binary separation around t — 0. There, the |i| 
term dominates, corresponding to the free parabolic so- 
lution derived by slightly varying the time at periapsis 



(eq. [T4"l) . The latter, however, does not diverge at t = 0. 
Instead, it describes a circle centered around x = with 
y = R, where R is the radius of the circle. 

A key observation for our problem is that all parabolic 
orbits which reduce to the A s (B s ) solutions are asym- 
metric (symmetric) function of time where the subscript 
s stand for x, y or z. We conclude that we can ferry a free 
solution across t = 0, from negative to positive times, by 
simply changing the sign of its As coefficients. 

A concern might arise from the fact that our assump- 
tion of small binary separation relative to the black hole 
may break down at times where the binary is close to 
its periapsis. Indeed, at small enough times, the binary 
separation is about a, | i | x / 3 while the distance to the 
black hole decreases as r m ~ a (M/m) 1 / 3 !*) 2 / 3 . They 
match at \t\ — (m/M) 1 / 3 , from which follows that, if 

the binary gets within a distance of r m /a < (m/M) 1 / 9 
from the black hole, our approximation is no longer valid. 
This occurs for a penetration factor for the orbit smaller 
than D < (m/M) 2 / 9 . When this happens, also the 
parabolic formalism incorrectly traces the trajectory for 
\t\ < (m/M) 1 / 3 . Specifically the true orbit deviates from 
the circle given by eq. (|14[) . Yet, the rest of our conclu- 
sions, including the energy of the particles will not be 
affected. Using the same work argument in our intro- 
duction, one derives that the energy gained during the 
periapse passage is smaller by a factor of D 2 than that 
gained around the tidal radius. One factor of D come 
from the smaller distance around the black hole, and the 
other from the larger distance between the stars resulting 
in a smaller force between them. We conclude that - for 
any impact parameter- both energy and orbit obtained 
with the method described above are fairly accurate for 
|t| > (m/M) 1 / 3 . 

4. NUMERICAL VERIFICATIONS 

We test our approximated equations against three- 
body simulations of a binary evolving around a much 
more massive black hole. 

Both the three-body code and the code that numeri- 
cally integrates our equations are provided with a 4th- 
order Runge-Kutta integration scheme. In the three- 
body code the binary center of mass is initially either 
on a parabolic ( §4. 1|) or on a radial ( §4.21) orbit. The 
binary consists of a primary star with mass mi and of 
a secondary with m,2 — mi/3. In the examples shown 
in this section, the black hole-binary mass ratio is set to 
be M/m = 10 6 . The binary's orbit is assumed to be ini- 
tially circular in the comoving frame of the binary center 
of mass. The initial configuration of the system is char- 
actered by three parameters. First, the initial distance of 
the binary center of mass to the black hole, rg- However, 
as long as a simulation starts at a large enough radius, 
fo ^ r t, the orbits are largely independent of it. In our 
runs we assume r$ ~ 10r t , which is sufficient for con- 
vergency. Second, the initial (at io = t(ro) < 0) binary 
phase, 0o- We parameterize it using the effective phase 
cf> at t = (i.e. at the periapsis passage for a parabolic 
orbit or at r = for a radial one), 

0o = S ujt + 4>, 
where u> is the (constant) angular velocity of the binary 



6 




>> o : 




Fig. 1. — The orbit of the secondary star in the primary comoving 
frame. The binary center of mass is initially on a parabolic trajec- 
tory. We assume a prograde binary orbit, M/m = 10 6 ,mi/rri2 = 
3, D = 0.1, <j> = 47r/5, and ro = 10rt. The three-body solution 
(solid line) is compared with our approximated solution (dashed 
line). Squared (three-body) and cross (our solution) marks show 
the position of the secondary at t = —3,0,3 and 15. Lengths and 
time are in units of (m/M) 1 / 3 r p and (rJ/GM) 1 ' 2 , respectively. 

at r > ft- Naturally, the actual phase at t = is in 
general different from cj>, due to of the black hole tidal 
forces. All angles are measured from the x-axis. Finally, 
we should specify the direction of rotation of the binary 
S, as viewed in the non-rotating frame: i.e. the relative 
orientation of the the angular momentum of the binary 
around the black hole and of a star around the binary 
center of mass. For a planar orbit there are two possibil- 
ities: the angular momenta are aligned, S = 1, in which 
case we call the orbit prograde, or they are anti-aligned, 
S = — 1, and we call it a retrograde orbit. In addition, for 
the parabolic case, the system has a fourth parameter: 
the penetration factor D. 

4.1. Parabolic Orbits 

In this section, we compare the evolution of the bi- 
nary stars obtained with the three-body code and our 
parabolic formalism (eqs. l8l fl"Tj) . We assume a prograde 
orbit for the binary, D — 0.1 and <j> — 4tt/5. The result 
is shown in Figure [U where we plot the orbit of the sec- 
ondary star in the comoving frame of the primary. After 
disruption around periapsis, the secondary star is cap- 
tured by the black hole on an elliptical orbit, while the 
primary is ejected from the system. Clearly, the three- 
body curve is accurately reproduced by our set of ap- 
proximated equations. The energy, for example, differs 
at a 0.1% level after the binary disruption. 

In addition, we compared the performance of the two 
numerical methods, for M/m ^> 1. The calculation 
of our approximated solution is faster for two reasons. 
First, for forth order scheme, numerical convergence is 
achieved with a time-step about ~ (M/m) 1 / 12 times 
smaller. Second, we integrate seven equations, instead 
of eighteen, since we only follow the binary stars rela- 
tive distance. This moderate speed-up, combined with 
the fact that most dependencies are analytic, allows us 
to explore more easily a wider portion of the parameter 
space. 



Fig. 2. — As Figure [T] but for an initial radial orbit for the binary 
center of mass. All relevant parameters are the same. Here, the 
position of the secondary is marked at t = — 10~ 2 , — 10 _1 , —1 and 
—3. Length and time are in unit of a and ^/a 3 /Gm, respectively. 

4.2. Radial Orbits 

We now consider a binary moving on a radial orbit 
of zero energy and test the accuracy of our radial ap- 
proximation (eqs. !26H28p . The trajectory of the separa- 
tion between the stars is shown in Figure for a bi- 
nary with a phase cf> — 4tt/5. Here, we plot the or- 
bit up to t = —2 x 10~ 3 . The assumption of a rela- 
tive small star separation is expected to become invalid 
around \t\ ~ (m/M) 1 / 3 = 10~ 2 . Again, there is a good 
agreement between our solution and the three-body cal- 
culation. Comparing the two terms on the right-hand 
side of eq. (|26|) . one expects that the tidal force domi- 
nates for \t\ <C 1. Indeed, we find numerically that the 
deviation from the initially circular orbit becomes signif- 
icant around \t\ ~ 1. For |t| <C 1, the orbit approaches 
the free solutions x cx |t| -1 / 3 and y cx It] 1 / 3 . Finally, we 
also reproduce at a percentage level the energy evolution 
as a function of time. For the specific example in Fig- 
ure [H we get an energy that differs only by half a percent 
at t ~ from the three-body result. 

4.3. Deep Penetrators 

In the previous sections, we have established that both 
our sets of equations well reproduce the full three-body 
calculation. We now investigate their relationship. In 
particular, we show how a parabolic solution reduces to 
a radial one, in the limit in which a binary plunges deeply 
into the tidal sphere of influence of the black hole. At 
this aim, we employ our parabolic formalism to numer- 
ically calculate binary orbits with increasingly smaller 
penetration factors, (Figure [3j solid lines). In this exam- 
ple, the orbit degenerate to a radial one (dashed line) for 
D ~ 10" 5 . 

5. RESULTS 

We are now in the position to statistically describe the 
properties of a binary after its encounter with a black 
hole. We focus here on the planar case. For given masses 
of the three-bodies and semi-major-axis of the binary 
stars, the fate and the final energy of the binary members 
depend only on the penetration factor D — r p /r t and on 
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Fig. 3. — As Figure [T] but for different penetration factors as 
labelled. The binary's orbit is prograde with parameters M/m = 
10 6 ,</!> = 0, ro = 10 tv Lengths are in unit of the initial binary 
separation. The parabolic approximation (solid curves) and the 
three-body results (dotted curves) are pratically indistinguishable. 
In addition, we plot the radial orbit (dashed curve) for the same 
parameters. 

the effective binary phase (j). 

5.1. The Fate of the binary 

Contrary to naive expectations, we find that for D <C 1 
a non-negligible fraction of the binaries are not disrupted 
(Fig. 2]). For D < 10 -1 , the fraction of disrupted cases 
saturates at a level of ~ 80%. The best chance of disrup- 
tion is for a binary in prograde orbit, with D = 0.15. A 
comparison between the curves for prograde (labelled P p ) 
and retrograde (P T ) orbits underlines the well known fact 
that retrograde binaries tend to be more stable against 
tidal disruption. Nevertheless, the two curves converges 
for D <C 1. This is because in this limit the binary cen- 
ter of mass approaches the black hole in an almost radial 
fashion, so that the angular momentum of the binary 
around the black hole is very close to zero. Indeed, in 
the radial formalism, there is no distinction between pro- 
grade and retrograde orbits. Consistently, when we cal- 
culate the disruption probability for a binary on a radial 
orbit, we find a survival fraction of 19%. Interestingly, 
we observe that binaries that avoid disruption tend to 
tighten, with their final semi-major axis af < a. For 
instance, for radial orbits, the encounter with the black 
hole produces harder binaries in about 90% of the cases. 

5.2. Final energy of ejected stars 

In Figure we show the contour plots for the final 
energy, at t 1 of a member of a binary in parabolic 
orbit that faces either capture by the black hole (E = 
constant < 0) or ejection (E — constant > 0) from the 
three-body system. The flat amaranth region is where 
binaries are not disrupted. For a binary approching the 
black hole on a prograde orbit (upper panel), a rather 
shallow penetration factor of D = 2.1 is already sufficient 
to be torn apart by tidal forces, while no disruption oc- 
curs for a retrogarde orbit (lower panel) with D > 0.44. 
We chose to plot the lower panel with reversed axises to 
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Fig. 4. — Probability of disruption as a function of the penetra- 
tion factor D, for prograde orbit (P p , red solid line), retrograde 
orbit (P r , blue solid line) and total P = (P p + P r )/2. The binary 
phase [0 n] is sampled with 3000 equally spaced grid points for 
each penetration factor. 

emphasize how they tend towards the same energy dis- 
tribution for D — > 0. The limiting distribution depends 
only on (f> and it is given, of course, by the energy plot 
obtained for radial orbits (Fig. [6]). Note again that in the 
radial limit there is a finite region (— 0.317T < <f> < — 0.127T 
and 0.69-7T < <f> < 0.887r) where binaries survive. 

Let's now consider the pattern of the energy contours 
(Figures [5] and [6|). For the same penetration factor, it 
shows how a change of ir in the effective phase translates 
into an energy which is equal in absolute value but oppo- 
site in sign. Therefore, as noted before, a uniform distri- 
bution in 4> implies that, when the binary is disrupted, 
each body has 50% chance of being ejected. Another 
noticeable feature is the presence and position of steep 
narrow peaks next to large plateaux. For prograde orbits, 
peaks of \E\ ~ 27 are attained not in the radial limit but 
for a finite penetration factor D ~ 10 _1 (see also solid 
red curve in the lower panel of Fig. [7|). On the other 
hand, for retrograde orbits, peaks -or rather "hills"- of 
a more modest maximum energy of \E\ ~ 1.8 emerge 
gradually as D gets smaller (see red dashed curve, same 
panel). Beside the maximal ejection energy, we also plot 
in Fig. [7] the ejection energy averaged over the binary 
phase, as a function of D (upper panel). In both plots, a 
rather large peak is present only for prograde orbits (red 
solid curves), and it is situated between 10 _1 < D < 1. 

A brief summary with the quantitative results of this 
section follows. 

• The largest D for which there is disruption is D = 
2.1 for prograde orbits and D — 0.44 for retrograde 
orbits. 

• The maximum energy is 27.3 for prograde orbits 
and 1.8 for retrograde orbits. 

• The average energy is 1.46 for prograde orbits, 0.50 
for retrograde orbits, 1.36 for prograde and retro- 
grade orbits together. In averaging, we have as- 
sumed that D, or equivalently r p , is uniformly dis- 
tributed. 




Fig. 5. — Final energy contour plots, in the phase-penetration 
factor plane for a given initial separation a. upper panel: pro- 
grade orbits, lower panel retrograde orbits. Energy is in unit of 
(Gmirri2/a)(M/m) 1 / 3 and the phase in unit of it. In the lower 
panel axis are reversed. The energies have been computed using 
the restricted parabolic formalism. The figure is constructed with 
a grid of 400 <j> values (equally spaced linearly) and 100 D values 
(equally spaced logarithmically) . We only plot the final energy of a 
member of a disrupted binary, which can be either captured by the 
black hole E < 0, or ejected E > 0. The flat amaranth (dark red) 
region shows the region in the parameter space where the binary 
survives disruption. 
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Fig. 7. — top panel: Ejection energy averaged over binary phase 
as a function of the penetration factor. The solid (dashed) line is 
for prograde (retrograde) orbits. The energies for a binary with a 
given a are plotted as red lines, in unit of (Gmim,2/a)(M/m) 1 ' 3 . 
The energies binaries with a given iJ m ; n are plotted as blue lines, in 
unit of (Gmim,2/-R m i n )(M/m) 1 / 3 . bottom panel: As upper panel, 
but for the maximum ejection energy for a given D, as a function 
of D. The average and maximum values are evaluated for the 
absolute value of the energy \E\. 

• The highest chance for disruption is for prograde 
orbits with D = 0.15. 




<|>[fl] 

Fig. 6. — Final energy as a function of for a member of a 
deep penetrating disrupted binary, with a given initial separation 
a. Energy is in unit of (Gmi 7712 /a) (M/m) 1 / 3 and the phase in unit 
of 7r. The red solid line is for a binary in radial orbit. Results for the 
parabolic formalism are shown for prograde orbits with D = 10 — 3 
(blue dotted), 10 -2 (green dashed-dotted) 3xl0 -2 (black dashed). 
The boundaries of the intervals in <j> where disruption of the binary 
occurs are marked with circles. The phase range is sampled with 
6000 equally spaced grid points. 



6. FINITE STELLAR SIZES 

Up until now, we have ignored the finite size of the 
stars. However, the highest energies or ejection velocities 
are obtained when stars, under the action of the black 
hole tides, get closer to each other before the binary dis- 
solves. Specifically, note the sharp peaks in Fig. [7] (lower 
panel) with exceptional high final energy. Inspection of 
these orbits, reveals that the minimal separation that 
these binaries attain over their evolution is more than 
ten times smaller than their initial separation. Those or- 
bits, however, are not always physical: if two stars start 
as almost contact binaries and contract further, they may 
collide and merge or tidally disrupt each other. We there- 
fore consider in the following a given minimum separation 
-Rmin ~~ about the sum of the radii of the two stars- below 
which a binary cannot shrink. 

The parabolic and radial formalisms discussed above 
provide the relative trajectory of the two stars, but the 
dimensions are arbitrary and can be rescaled. For each of 
these trajectories, we can find the minimum dimension- 
less distance ^/ x 2 + y 2 + z 2 | m in over the whole orbit, and 
scale it to be equal to i? m i n . This amounts to deriving 
the initial binary separation results in an orbits whose 
minimal separation is i? m i n : 

a = R m i n /D^/ x 2 + y 2 + z 2 | m in in the parabolic for- 
malism, and a = R m - ln / \J x 2 + y 2 + z 2 | m i n in the radial 
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Fig. 8. — As Fig.[5]but for a given minimal distance i? m i n . Here, 
energy is in unit of (Gmim2 / fl m ; n ) (M/m) 1/13 and the phase in 
unit of 7T. Note the energy scale, from —2 to 2, a factor of 5 
smaller than the energy scale needed for prograde orbits with a 
fixed a (Figure [5] upper panel). 

formalism. 

Assuming this a as our initial condition, the energy is 
given by 



Ei — — E2 



Gm\m,2 ( M 



1/3 



— y/x 2 + y 2 + z 2 \, 
m 



(1 + cos/) 5 



(xcos/ + y sin/) + 



(35) 

— sin fx + (1 + cos f)y 
V2 



or by 



1/3 



E x = -E 2 = — ) y / x 2 +y 2 + z 2 \ 1 



2/3 



(36) 



|t|~ 4 / 3 (s + 3tx) , 



in our parabolic and radial formalism, respectively. 

Our analytic approach, therefore, allows us to recast 
our results in terms of a given minimal distance between 
the stars, (Figs. [5] and [5]), rather than a given initial 
separation, (Figs. [5] and [S]). In comparison, the main 
notable feature is that the scale of Fig. [5] is much nar- 
rower, i.e., the energy is more sensitive to i? m i n than to 
a. This is shown more clearly in Fig. where the mean 
and especially the maximal energy for a fixed i? m in are 
flatter function of D. The maximal energies in units of 
(Gm 1 m 2 /]? mm ) (M/m) 1 / 3 are: 

• 1.64 for prograde orbits (obtained for D = 0.12) 
and 1.60 for retrograde orbits (obtained for D = 0). 
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Fig. 10. — The ejection velocity averaged uniformely over phase 
as a function of the penetration factor. The average is over 3000 
<f> equally spaced, upper panel: The ejection velocity for a given 
a in unit of yj2Gm 2 /a{M/m) 1 / & . We plot it for prograde orbits 
only (red solid line, (V) ), for retrograde orbits only (blue solid 
line, (V) ) and the weighted average over these two orbital direc- 
tions (black dashed line, (V) T ). See text for details. Lower panel: 
as upper panel, but for a given minimum possible distance R m i n 
between stars: here (V) is in unit of 2Gm2 / Rmin (M/m) 1 /® . 

7. EJECTION VELOCITIES 

In the quest for HVSs the observed quantities are the 
radial velocity and the star mass. Neglecting the Galac- 
tic potential, an ejected star (E\ > 0) with mass mi 
has a velocity at infinity of v\ = yl 1E\ J m\ ■ A compre- 
hensive study of velocity and mass distributions, under 
various assumptions for the binary population is deferred 
to a subsequent paper. Here, we calculate the mean ve- 
locity for a given binary, assuming that the binary plane 
is the same as the orbital plane around the black hole 
and that the binary phase is uniformly distributed. The 
results are shown in Fig. [10] (upper panel), as a fucntion 
of the penetration factor. If we call (v) , the average ve- 
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locity for a prograde orbit (red solid line) and (v) r that 
for a retrograde orbit (blue solid line), then the total 
average is (u) tot = n p (v) p + n r (u) r (black-dashed line), 
where n r and n p are the fraction of ejected stars that 
originate from retrograde and prograde binaries. Clearly 
there is no gain in deeper impacts, where the mean sat- 
urates at ~ 0M-J2Gm2~/a(M /m) 1/6 . Higher velocity 
-but only by 30%- can be attained closer to the tidal 
radius. We also calculate similar means, but for bina- 
ries that, under the black hole tides, shrink to the min- 
imum relative distance i? m i n (see $6}. Harder binaries 
- which would be responsible for the higher mean ve- 
locities - would now coalesce instead. This results in an 
even flatter mean velocity as a function of D, which tends 
to - 0.94 v /2Gm 2 /i? mi „ (M/m) 1 / 6 for D -> 0, (Fig. M 
lower panel). 

8. DISCUSSION AND CONCLUSIONS 

The ultimate goal of our work is to statistically car- 
acterize the population of stars originated from tidally 
disrupted binaries, and compare it with observations of 
hyper velocity stars. 

To this purpose, we derive in this paper the equations 
of motion and energy for a member of a close binary, 
which suffers an encounter with a third far more massive 
body. For the general case, we assume a parabolic orbit 
for an effective center of mass of the binary (Sj2|. This 
is in contrast with previous works, (e.g. iBromlev et al.1 
2006) who considered binaries on hyperbolic or ellipti- 
cal orbits. However, our assumption is justified since the 
orbits of binaries that are candidates to produce hyper 
velocity stars are very eccentric. The periapsis passage 
of a very eccentric orbit could be modeled by a parabolic 
orbit to a good a pproximation. This is in agreement 
with the finding of Bromley et al.l (|2006l ) that the initial 
binary velocity is of little influence on the final outcome. 
Nevertheless, a binary on a hyperbolic orbit has a total 
positive energy, and allows for a new disruption channel, 
in which both stars are ejected. This channel is rela- 
tively rare, since the energy of the binary is typically 
small compared to the typical ejection energies. More- 
over, such double ejections could not lead to hyperve- 
locity stars, since the energy is limited to the original- 
small- positive energy of the binary. Our formalism with 
zero total energy (except the small negative binary bind- 
ing energy), does not allow for double ejections. 

A parabolic trajectory for the binary approaches a ra- 
dial one for a very close encounter with the massive body. 
This observation leads us to explicitly adopt a radial or- 
bit with zero energy, in order to follow the limiting case 
of a deep penetrator (Q. This simpler set of equations 
allow us to easily trace a close encounter, that otherwise 
would require high accuracy when calculated with a full 
three-body code. 

Our formalism can be applied quite generally to ex- 
plore the fate of a binary with arbitrary orbital parame- 
ters. However, in this paper, we only focused on results 
for circular coplanar binaries. The inclination of the bi- 
nary, as well as its eccentricity, is expected to affect our 
results quantitatively. We reserve such a study to a forth- 
coming paper. Nevertheless, we can already reach some 
conclusions and note quantitative differences with results 
found in literature. 



The main feature of HVSs is of course their unusu- 
ally high radial velocity. For equal mass stars, mi = 
rri2 = mi t 2, the expression Gmiiri2 / R m in is simply 
"ti^Vesc/^ where v esc is the escape velocity from the 
surface of the stars. Since the maximal ejection en- 
ergy is 1.6 Gm\m2/R m in{M/m) 1 ^ (Fig. [7]), we derive 
for the Galactic black hole a corresponding velocity of 
0.9 Ves^M/m) 1 / 6 = 9t> csc . The escape velocity for main 
sequence stars is about 600-800 km/s in the mass range 
of 1 — 10M Q . Therefore velocities can be as large as 
~ 7000 km/s even for binaries limited to the plane. Scal- 
ing f or the different masses (and thus iw ) assumed by 
iHillsl (fl988l ) and bv IBromlev eTaTl (pool , we find that 
a velocity of 0.9 v esc (M/m) 1 ^ 6 is respectively a factor of 
~ 1.2 and ~ 1.8 higher than their maximum value of 4000 
Km/s. Even higher velocities can be acheived when the 
binary mass ratio is large and the lighter star is ejected. 
Then, the ejected star travels with maximal velocity of 
1.3 Ves^M/m) 1 / 6 , which is around ~ 10,000 km/s for 
our Galactic center black hole where M/m ~ 10 6 . 

For a comparison with observations, it is also impor- 
tant to determine the ejection probability. IBromlev et al.l 
(2006) find that the probability for disruption as a func- 
tion of the penetration factor goes roughly as P c) ~ 
1 — D/2.2 . The maximum D for disruption is there- 
fore very similar to ours, D max ~ 2.1 (see Fig. [10]) . In 
contrast, we remark a quantitatively different behaviour 
of the probability function for binaries that approach 
closely the black hole, D — ► 0. Specifically, their P C j — > 1 
implies that all binaries in this limit are disrupted. This 
may look intuitively sound since binaries that penetrate 
so deeply experience very strong tides. However, as we 
stressed in SJH the stars in such binaries separate, but 
approach each other again on the way back from the 
black hole. In the planar case we find that ~ 20% of 
the deeply penetrating binaries survive (Fig. g]). Even 
when taking a uniform distribution for the inclination 
of the binary plane into account, the percentage of sur- 
vival is still 10% for close encounters. In the planar case, 
Figure H] indicates that the disruption probability is, in 
fact, not monotonic in D: for prograde orbits, it peaks 
at D = 0.15 and it is almost unity (98.0 - 98.8%) for 
0.06 < D < 0.3. 

For large mass ratios mi/mz > 10, IBromlev et alJ 
(2006) find that the heavier has consistently more chance 
to be ejected. We showed here that the probability is 
50%. This result is not limited to zero inclination of ec- 
centricity. This fact together with the rarity of a very 
massive star somewhat weakens the claim that the star 
SO-2 was created by a disruptio n of a binary, in which a 
60 M & companion was ejected (jGould fc Quillenl feOOS). 
This conclusion was based on the observed orbital pa- 
rameters of SO-2. However, the timescale to significantly 
change its shor periapsis distance may quite short, and 
a careful study of the dynamical processes in the inner 
regions of the Galactic Center is needed to asses it. 
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